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As a first step towards understanding the morphology of PdO crystals we performed a systematic 
full-potential density-functional theory study of all possible (1 x 1) terminations of the low-index 
surfaces of tetragonal PdO. Applying the concept of first-principles atomistic thermodynamics we 
analyze the composition, structure and stability of these PdO orientations in equilibrium with an 
arbitrary oxygen environment. Within the studied subset of (1 x 1) geometries the polar PdO- 
terminated PdO (100) orientation turns out to be surprisingly stable over the whole range of exper- 
imentally accessible gas phase conditions. Setting up a constrained Wulff construction within the 
compiled data set, this PdO(100)-PdO facet correspondingly dominates the obtained polyhedron by 
far. The real PdO crystallite shape will however likely be affected by surface reconstructions, which 
are not covered by the present study. 



I. INTRODUCTION 



Metal oxides are compounds of widespread technolog- 
ical interest, and one field of application is catalysis, 
where they can act as the active material or the (often 
not that passive) support 1 '?. To obtain a microscopic un- 
derstanding of the function of these compounds in such 
applications it is necessary to know their surface atomic 
structure, which is also influenced by temperature and 
partial pressures in the environment. This can be par- 
ticularly important for oxygen containing environments, 
where the stability of different surface terminations of 
varying stoichiometry may well be anticipated as a func- 
tion of oxygen in the surrounding gas phase. Consider- 
ing the technological importance of oxides the scarcity 
of such atomic-level information even for well-ordered 
single-crystal surfaces is surprising. The little that is 
known stems almost exclusively from ultra-high vacuum 
(UHV) experiments, and is furthermore largely concen- 
trated on some specific oxides like the vanadium oxides, 
rutile (Ti02) or corundum (A^O^)^. For other oxides 
often not even the low-energy surface orientations are 
firmly established, and this also applies for the case of 
palladium oxide (PdO). 

Although PdO is renowned for its high activity in 
the catalytic combustion of methana^^, and the in- 
volvement of oxidic structures in high-pressure CO oxi- 
dation reactions at Pd surfaces is now being discussed^, 
virtually no information about the electronic and ge- 
ometric structure of PdO surfaces is presently avail- 
able, neither from experiment nor from theory. As a 
first step we therefore investigate the surface structure 
and composition of all low-index surfaces of tetrago- 
nal PdO in equilibrium with an arbitrary oxygen envi- 
ronment using the concept of first-principles atomistic 
thermodynamic&l&H'H based on density- functional the- 
ory (DFT) calculations (Section II). Lacking any experi- 
mental data on surface reconstructions we first focus on 
all possible (lxl) terminations and set up a constrained 
Wulff construction for this limited set for the whole range 



of experimentally accessible gas phase conditions (Sec- 
tion IIIB+C). This provides a first data base against 
which future models of reconstructed surfaces may be 
compared, in particular by how much they would have 
to reduce the surface free energy in order to have corre- 
sponding facets contribute significantly to the real PdO 
crystal shape. In addition, relatively high surface free 
energies and work functions obtained for (lxl) termi- 
nations might point at likely candidates for surface re- 
constructions. Interestingly, within the studied subset 
of (1 x 1) geometries, one termination (PdO(lOO)-PdO) 
turns out to be much more stable than all others, and 
correspondingly dominates our constrained Wulff poly- 
hedron by far. And this although it represents a so-called 
polar termination, which are traditionally dismissed on 
electrostatic groundsi&i4 (Section HID). 



II. THEORY 

A. Atomistic thermodynamics 

In order to describe the thermodynamic stability of 
PdO surfaces in an oxygen environment we use DFT 
total-energy calculations as input to atomistic thermody- 
namics considerations2ii2iiiii2ii£ii&, which treat the effect 
of the surrounding gas phase via the contact with a corre- 
sponding reservoir. In equilibrium with such a reservoir 
the most stable surface structure in the constant temper- 
ature and pressure (T,p)-ensemble minimizes the surface 
free energy, which is defined as 



(1) 



Here G is the Gibbs free energy of the solid with the sur- 
face of area A, to which in a supercell calculation with 
symmetric slabs both the top and bottom surface con- 
tribute equally. in(T,pi) is the chemical potential of the 
various species i present in the system, i.e. in this case 
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i = Pd, O. Ni gives the number of atoms of the ith com- 
ponent in the solid. For not too low temperatures and 
sufficiently large particles bulk PdO may further be as- 
sumed as a second thermodynamic reservoir with which 
the surface is equilibrated. This constrains the chemical 
potentials of O and Pd to the Gibbs free energy of bulk 
PdO, ffpjo (where the small g denotes the Gibbs free 
energy per formula unit), and allows to eliminate /xpd 
from eq. . The remaining quantities to be determined 
for the calculation of the surface free energy are then the 
chemical potential of the oxygen gas phase, as well as the 
difference in Gibbs free energies of slab and bulk PdO. 

The computation of the prior is straight-forward, as 
fio is completely fixed by the surrounding gas phase 
reservoir, which may well be approximated as an ideal 
gas. Ideal gas laws then relate the chemical potential 
to temperature and pressur o 15 i 17 , and we will convert 
the dependence of the surface free energy on fio (T, p) in 
all figures also into the more intuitive pressure scales at 
T = 300 K and T = 600 K. The second input to j(T,p), 
i.e. the Gibbs free energy difference of the bulk phase 
and the slab, receives contributions from changes in the 
vibrational and configurational degrees of freedom at the 
surface, as well as from the pV-teim and as leading con- 
tribution from the difference in total energies. From a di- 
mensional analysis, the pV-term. may well be neglected^. 
The configurational contribution for a system like PdO 
can further be estimated as below 5meV/A 2 ji& again 
negligible for a study that aims at a first, rather coarse 
comparison of different (lxl) surface terminations. 

The vibrational contribution can be obtained from 
first-principles using the computed phonon density of 
states (DOS) at the surface and in the bulk, cf. e.g. 
ref. Alternatively, the Einstein approximation to the 
phonon DOS can be employed to get an order of magni- 
tude estimate^. Allowing a 50% variation of the charac- 
teristic frequency^ of each atom type at the surface, the 
vibrational contribution to the surface free energy at all 
considered PdO terminations stays in this model always 
within a range of about 10-20 meV/A 2 for all tempera- 
tures up to T = 600 K. Although this is certainly not a 
small contribution in general anymore, we will neverthe- 
less neglect it in this particular study. Being interested 
in a coarse, first screening of the stability of various PdO 
terminations, it will become apparent below that even 
a 10-20 meV/ A 2 contribution will not affect the physical 
conclusions drawn. If in other studies a higher accuracy 
is required, the vibrational term may either be taken into 
account by e.g. simplified treatments of the most rele- 
vant vibrational modest, or eventually by a full DFT 
calculation of the phonon DOS. Here, we will neglect 
it however, and may with all other approximations dis- 
cussed in this section then replace the difference of bulk 
and slab Gibbs free energies entering into the computa- 
tion of j(T,p) simply by the corresponding difference of 
total energies. 



B. DFT basis set and convergence 

The DFT total energies that are thus needed are ob- 
tained using a mixture of the full-potential augmented 
plane wave + local orbitals (APW+lo) and the linear 
augmented plane wave (LAPW) method together with 
the generalized gradient approximation (GGA) for the 
exchange-correlation functional^! as implemented in the 
WIEN2k codo 22 i 23 i 24 . To simulate the different PdO sur- 
faces we use supercells containing symmetric slabs with 7- 
11 layers and 12-15 A vacuum between subsequent slabs. 
The outermost 2-4 layers are fully relaxed for all surfaces. 
For all orientations we have also performed test calcula- 
tions with thicker slabs and relaxing more layers, without 
obtaining any significant changes to the chosen setup (< 
3meV/A 2 in the surface energy). 

The parameters for the mixed APW+lo and LAPW 
(L/APW+lo) basis set are: R^ T = 1.8bohr, R^ T = 
1.3bohr, wave function expansion inside the muffin tins 
up to l^ ax = 12, potential expansion up to l^ x = 6. 
The Brillouin zone (BZ) integration is performed using 
Monkhorst-Pack (MP) grids with 6, 8 and 16 k-points 
in the irreducible part of the BZ for the (110) / (111), 
the (100)/(010)/(101)/(011) and the (001) directions, re- 
spectively. The energy cutoff for the plane wave repre- 
sentation in the interstitial region is -E^lx = 17 Ry for 
the wave function and E^ x = 196 Ry for the potential. 
With these basis sets the surface energies of the different 
PdO surfaces are converged within 1-2 me v/A 2 regarding 
the k-points and 3-4meV/A 2 regarding E^ ax . Errors in- 
troduced by the fundamental approximation in the DFT 
approach, i.e. our use of a GGA as exchange-correlation 
functional, will be commented on below. 



III. RESULTS 

Since there is virtually no atomically-resolved informa- 
tion about the structure and composition of crystalline 
PdO surfaces available from the experimental side, we 
start our theoretical investigation from a very basic point 
of view. That is, to get a first idea about the geometric 
and electronic properties of different PdO surfaces we fo- 
cus here on a rather coarse comparison of the subset of 
all possible (lxl) terminations of the low- index surfaces 
of tetragonal PdO. 

A. Geometric bulk and surface structure 

PdO crystallizes in a tetragonal structure within the 
space group -Df h — • There are two formula units of PdO 
in the tetragonal unit cell with Pd atoms at all corners 
and in the centre, and O atoms at (0, 1/2, 1/4), (0, 1/2, 
3/4), (1, 1/2, 1/4) and (1, 1/2, 3/4), as shown in Fig. 
n All Pd resp. O atoms are equivalent, with each Pd 
atom planar coordinated by 4 oxygen atoms, and each O 
atom tetrahedrally surrounded by 4 Pd atoms. Within 
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FIG. 1: The tetragonal bulk unit cell of PdO. Small dark 
spheres indicate oxygen atoms, large light ones Pd atoms. 



our DFT-GGA approach the optimized lattice constants 
of the PdO unit cell are obtained as a = 3.051 A and 
c = 5.495 A, which is in reasonably good agreement with 
the experimental lattice constants, a oxp = 3.043 A and 
'-•cxp — 5.336 J^^^". 

Due to the tetragonal structure of the PdO unit cell, 
there are 5 inequivalent low-index orientations, each with 
2-3 different (1 x 1) surface terminations. The PdO(lOO) 
surface (parallel to the yz-plane) is equivalent to the 
(010) surface (parallel to the xz-plane). For this surface 
direction there are two different (lxl) surface termina- 
tions, one containing only Pd atoms in the topmost layer 
(PdO(lOO)-Pd), and the other Pd as well as O atoms 
(PdO(lOO)-PdO). Schematic drawings of these and all 
following surface geometries are shown as insets in Fig. 
121 below. The PdO(OOl) surface is parallel to the xy- 
plane, and there are again two (1 x 1) terminations, one 
with only Pd atoms (PdO(OOl)-Pd) and one with only O 
atoms (PdO(OOl)-O) in the topmost layer. 

For the PdO(lOl) (equivalent to PdO(Oll)) surface 
there exist three different terminations. One termina- 
tion is stoichiometric cutting the stacking of 0-Pd2-0 
trilayers along the (101) direction just between consec- 
utive trilayers (PdO(lOl)). The other two terminate ei- 
ther after the Pd layer (the O-deficient PdO(lOl)-Pd) or 
have two O layers at the top (the O-rich PdO(lOl)-O). 
The remaining PdO(llO) and (111) directions are on the 
other hand characterized by alternating layers of Pd and 
O atoms along the surface normal, exhibiting therefore 
either a completely Pd (PdO(llO)-Pd, resp. PdO(lll)- 
Pd) or a completely O (PdO(llO)-O, resp. PdO(lll)-O) 
terminated surface. 

Looking at these 11 different (lxl) terminations we 
notice that only one of them is stoichiometric. The 
other 10 exhibit either an excess of oxygen or palladium 
atoms, and belong thus to the class of so-called polar 
surfacesiii4, the stability issue of which will be discussed 
in more detail in Section HID. 



B. Surface free energies 

As mentioned above we want to analyze the stability 
of these different PdO surfaces when in contact with an 
oxygen gas phase characterized by a given O chemical po- 
tential. This [io(T,p) can experimentally (and assuming 
that thermodynamic equilibrium applies) only be varied 
within certain boundaries. The lower boundary, which 
will be called the O-poor limit, is defined by the decom- 
position of the oxide into palladium metal and oxygen. 
A reasonable upper boundary for /zo on the other hand 
( O-rich limit) is given by gas phase conditions that are so 
oxygen-rich, that O condensation will start on the sam- 
ple at low enough temperatures. Appropriate and well- 
defined estimates for these limits are given by 1 - 



AG f (T = 0K, P - 0) < A^ (T,po 2 ) < 



(2) 



where the O chemical potential is referenced with re- 
spect to the total energy of an oxygen molecule, A/xo = 
HO - (l/2)£^° tal , and AG/ (T = 0K,p = 0) is the low 
temperature limit for the heat of formation of PdO. 
Within DFT-GGA we compute -0.88 eV for this quan- 
tity, which compares well with the experimental value of 
AGf p (T -> 0K,p = latm) = -0.97 eV^. To also con- 
sider the uncertainty in these theoretically well-defined, 
but approximate limits for Afio, we will always plot the 
dependence of the surface free energies also some tenths 
of eV outside these boundaries. From this it will become 
apparent below that the uncertainty in the boundaries 
does not affect at all our physical conclusions drawn. 

We proceed to show all surface free energy plots of the 
eleven discussed (1 x 1) terminations in Fig. [2] Termi- 
nations with an O excess (deficiency) show a negative 
(positive) slope, i.e. they will become the more favorable 
(unfavorable) the more O-rich the surrounding gas phase 
is. Comparing the five plots shown in Fig. [2] the very 
high stability of the PdO(100)-PdO termination is eas- 
ily recognized. The only other terminations that exhibit 
comparably low surface energies are the PdO(lOl), and 
towards the O-rich limit the PdO(110)-O and PdO(lll)- 
O surfaces. All other considered terminations are rather 
high in energy, in particular all those that feature Pd 
atoms in their outermost layer. Looking at the energy 
scale in the plots we also notice that this stability or- 
dering will not be affected by the afore discussed 10- 
20meV/A 2 uncertainty in our surface free energies. 

As already indicated above this uncertainty in the 
surface energies does not yet include the error due to 
the choice of GGA as exchange-correlation functional. 
We have therefore calculated the surface free energies 
of all considered terminations also within the local den- 
sity approximation (LDA)21. To set up the supercells 
for these LDA calculations we first optimized the lat- 
tice constants for PdO bulk, obtaining olda = 2.990 A 
and clda = 5.292 A, i.e. values that are as expected 
slightly smaller than the afore mentioned experimental 
lattice constants. After a full relaxation of the 2-4 out- 
ermost layers, the relative changes in the surface geome- 
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FIG. 2: Surface free energies for the 5 inequivalent low-index PdO surfaces. Solid lines indicate (1 x 1) terminations, and dashed 
lines larger unit cell reconstructions discussed in Section HID. The vertical dotted lines specify the range of A/io considered in 
this study (see text), while in the top two x axes the dependence on the O chemical potential has been converted into pressure 
scales at T — 300 K and T — 600 K. The insets show the surface geometries of the corresponding (lxl) terminations (O small 
dark spheres, Pd large light spheres). 



TABLE I: Surface free energies of all low-index PdO (1 x 1) 
terminations at the oxygen-poor limit, as calculated within 
the GGA or the LDA. All energies are in meV/A 2 , and the 
numbers in brackets denote the energetic difference with re- 
spect to the lowest-energy PdO(100)-PdO termination. 



Surface 
termination 




TO-poor 

GGA 




TO-poor 

LDA 


PdO(100)-PdO 


33 


(0) 


59 


(0) 


PdO(100)-Pd 


119 


(+86) 


170 


(+1H) 


PdOfOOlVO 


119 


f+86") 


162 


(+103) 


PdO(001)-Pd 


156 


(+123) 


212 


(+153) 


PdO(lOl) 


57 


(+24) 


86 


(+27) 


PdO(101)-O 


134 


(+101) 


180 


(+121) 


PdO(101)-Pd 


128 


(+95) 


173 


(+114) 


PdO(110)-O 


86 


(+53) 


119 


(+60) 


PdO(110)-Pd 


137 


(+104) 


173 


(+114) 


PdO(lll)-0 


72 


(+39) 


109 


(+50) 


PdO(lll)-Pd 


105 


(+72) 


143 


(+84) 



tries with respect to these bulk values are found to be 
very similar to the ones obtained within the GGA. The 
resulting absolute surface free energies of the eleven ter- 
minations in the oxygen-poor limit are listed in Table 
U (together with the corresponding GGA values). Al- 
though the absolute values of the LDA surface energies 
are 30-50meV/A 2 higher than within the GGA, the rela- 
tive differences between them are within 10-30 me v/A 2 . 
This is exemplified in Table [i] by also indicating the rel- 
ative energetic difference with respect to the PdO(100)- 
PdO termination, which is the lowest-energy surface in 
both the LDA and the GGA. As only this energetic order- 
ing among the considered (lxl) terminations enters into 
the targeted constrained Wulff construction, and we find 
this ordering to be similar for calculations with two quite 
different exchange-correlation functionals, we expect the 
DFT accuracy to be rather high for this system. 



C. Constrained Wulff construction 

With the results obtained for the surface free energies 
of the different (1 x 1) PdO terminations we set up a con- 
strained Wulff construction 2 ^ for a PdO single crystal. 
Since this construction is constrained to reflect only the 
studied (1 x 1) terminations, its intention is more to com- 
pare the relative energies of different surface orientations, 
rather than to really predict the equilibrium PdO crys- 
tallite shape (which will most likely be affected by sur- 
face reconstructions on which we comment below) . Since 
the surface free energies depend on the oxygen chemical 
potential, the obtained construction will vary with the 
conditions in the surrounding gas phase. We therefore 
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FIG. 3: Constrained Wulff construction at the oxygen-poor 
(left) and oxygen-rich limit (right). The construction is con- 
strained to the studied (1 x 1) terminations and reflects there- 
fore only the relative energies of different PdO orientations, 
rather than the real PdO crystal shape, which will most likely 
be affected by surface reconstructions. The polyhedra are 
symmetric with respect to the xy-plane, and only the upper 
half is shown correspondingly. 



present in Fig. [3] the Wulff polyhedra for the two consid- 
ered extremes, i.e. in the O-poor and the O-rich limit. 
Due to the tetragonal symmetry of PdO, the polyhedra 
are symmetric with respect to the xy-plane indicated in 
Fig. ^ and only the upper half of each polyhedron is 
shown correspondingly. 

As one would already assume from the stability plots in 
Fig. the low -energy PdO(100)-PdO surface forms the 
largely dominating facet (rectangular, dark gray facets) 
at both the O-poor and the O-rich limit of A/io- The 
other triangular, light gray facets are built by the stoi- 
chiometric PdO(lOl) surface termination. All other in- 
vestigated terminations do not contribute at all to the 
present construction, since the corresponding planes lie 
outside of the polyhedra and do not cross them at any 
point. The polyhedra in Fig. [21 are scaled in a way, that 
the absolute value of the area belonging to the PdO (101) 
termination is the same in the O-poor and O-rich limit, 
since the surface free energy of this termination is also 
constant with respect to A^o- I n turn, the remaining 
surface area built up by the O-rich PdO(100)-PdO termi- 
nation increases strongly in going from O-poor to O-rich 
limit. In the prior limit, this termination forms already 
72% of the whole area of the polyhedron, while at the 
latter limit this fraction has increased to even 94%. As 
would already be expected from the similar relative en- 
ergy differences, these Wulff construction properties are 
also almost equally obtained within the LDA: again the 
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TABLE II: Minimum energy by which surface reconstructions 
at the various facets would have to lower the surface free en- 
ergy, in order for the facets to touch the presently obtained 
constrained Wulff polyhedron (touching). Additionally, the 
corresponding lowering required for the facet to cover approxi- 
mately 10% of the total surface area of the polyhedron is listed 
(10% area). All energies in meV/A 2 (and percent changes) 
are given with respect to the lowest-energy (1 x 1) termination 
of the corresponding orientation. 



Orientation Touching 10% area 



Oxygen-poor: 

(001) -54 (-45%) -67 (-56%) 

(110) -40 (-47%) -44 (-51%) 

(111) -7 (-10%) -16 (-22%) 

Oxygen-rich: 

(001) -7 (-10%) -42 (-58%) 

(110) -39 (-80%) -41 (-84%) 

(111) -1 (-2%) -20 (-37%) 



PdO(100)-PdO and PdO(lOl) terminations are the only 
ones contributing to the polyhedron, and the PdO(100)- 
PdO facet covers 66% (82%) of the surface area in the 
O-poor (O-rich) limit. Although such a comparison be- 
tween two functionals forms no formal proof, we would 
therefore assume the obtained shape to be quite accurate 
with respect to this DFT approximation. 

This does, however, not comprise the major limitation 
of our study, given by the restriction to (1 x 1) termina- 
tions. Surface reconstructions could considerably lower 
the surface free energy of any of the PdO surfaces and 
correspondingly significantly influence the overall shape 
of the Wulff polyhedron. Unfortunately, to the best of 
our knowledge no experimental information on such sur- 
face reconstructions is presently available for crystalline 
PdO. Without any such information, not even on the sur- 
face periodicity, the phase space of possible reconstruc- 
tions is simply too huge to be assertively screened by 
todays first-principles techniques alone^S. Until such in- 
formation becomes available from experiment, the best 
we can do at the moment, is to check by how much such 
surface reconstructions would have to lower the surface 
free energy in order to give rise to significant changes in 
the real PdO crystal shape compared to the constrained 
construction presented in this study. This check is done 
for the three orientations presently not contributing to 
the exterior of the polyhedron, i.e. the (001), (110) and 
(111) facets, and since the surface reconstructions are of 
unknown stoichiometry their effect could be different in 
the O-poor and in the O-rich limit. 

Table m lists the corresponding values for each facet, 
of how much the surface free energy would have to be 
lowered by a reconstruction (with respect to the lowest- 
energy (lxl) termination presently considered in our 
constrained Wulff construction), such that the facet just 
touches the current polyhedron. Additionally given is the 




PdO(110)-O (lxl) PdO(110)-O (2x1) 

FIG. 4: Unit cells for stoichiometric reconstructions of the 
PdO(001)-O and PdO(110)-O terminations, achieved by sim- 
ply removing the oxygen atoms marked with crosses. The two 
left figures show configurations, where the oxygen atoms are 
taken out along a row; the right ones, where the oxygen forms 
a checkerboard pattern. 



corresponding value, required to have the reconstruction 
really contribute significantly to the total Wulff shape, 
which we consider to be the case when the facet covers 
approximately 10% of the total surface area. From the 
compiled data in Table ^ it seems that rather massive 
reconstructions would be necessary (both in the O-poor 
and in the O-rich limit) to have appreciable (001) and 
(110) facets in the PdO equilibrium crystal shape. A 
much smaller energy reduction would on the other hand 
be required to have a (111) facet contribute. 



D. Stability of PdO surfaces 

As already mentioned above only one of the 11 possi- 
ble (lxl) terminations, namely the PdO (101), is sto- 
ichiometric, whereas all others are so-called polar sur- 
faces, which are not expected to be stable on electrostatic 
groundsiii4. For the three different surface terminations 
in the (101) direction this ionic model certainly com- 
plies with our results, since the stoichiometric PdO(lOl) 
termination turns out much more stable than the po- 
lar PdO(101)-O and PdO(101)-Pd terminations. For the 
other orientations it is on the other hand not possible at 
all to truncate the tetragonal PdO structure in a (1 x 1) 
cell and achieve charge neutrality; and one might wonder 
whether this is the reason why e.g. the considered (001) 
and (110) (1 x 1) terminations exhibit such high surface 
energies? As most obvious guess we therefore performed 
fully relaxed calculations in larger unit cells, where we 
simply removed half the O atoms in the top layer in order 
to achieve stoichiometric terminations. As shown in Fig. 
21 there are two possibilities for both the (001) and the 
(110) orientation to remove these oxygen atoms, either 
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TABLE III: Work functions for the different (1 x 1) PdO 
surface terminations in eV. 



O-terminated 


$ (eV) 


Pd-terminated 


$ (eV) 






PdO(100)-Pd 


4.0 


PdO(001)-O 


7.9 


PdO(001)-Pd 


4.8 


PdO(101)-O 


7.7 


PdO(101)-Pd 


4.5 


PdO(110)-O 


7.2 


PdO(110)-Pd 


4.4 


PdO(lll)-0 


5.9 


PdO(lll)-Pd 


4.7 


PdO-terminated 


$ (eV) 


stoichiometric 


$ (eV) 


PdO(100)-PdO 


6.4 


PdO(lOl) 


5.4 



w 
•- 



s 

■5 
s 

£ 



X x 

XL- 



3-fold- 
coordinated 




(100) (101) (111) (110) (001) 



all along a row or in a checkerboard pattern. The corre- 
sponding surface free energies are drawn as dashed lines 
in Fig. |2 and are not at all lower than the corresponding 
polar (1 x 1) O-rich terminations. These results are there- 
fore at variance with the suggestion made by Ciuparu et 
al, that such a simple removal of O atoms should lead 
to charge compensated and thus stable PdO(OOl) and 
PdO(llO) surfaced. 

Apparently, charge neutrality is not the dominant fea- 
ture determining the energetic stability of PdO surfaces, 
as is also directly reflected in the very low surface free 
energy of the polar PdO(100)-PdO termination. This 
points at the most obvious shortcoming of the electro- 
static model, namely the assumption that all atoms of 
one species are identical and in the same charge state, 
irrespective of whether they are at the surface or in the 
bulk. As we had shown befor o 15 ! 30 , structural and elec- 
tronic relaxation at the surface allows for appreciable de- 
viations from this picture, such that other factors (like 
an appropriate excess stoichiometry at O-rich conditions) 
might well overrule the polarity issue. 

Still, that there is a different degree of polarity associ- 
ated with the different terminations, is nicely reflected in 
the corresponding work functions, cf. Table ITTT1 While 
the work function of the stoichiometric PdO(lOl) ter- 
mination is with 5.4 eV in a medium range, the work 
functions of all Pd-terminated surfaces are about 1.0 - 
1.5 eV lower. The ones of the O-terminated surfaces are 
on the other hand about 1.0 - 2.5 eV higher, just as ex- 
pected from the ionic model. Even the comparably low 
work function of the O-terminated PdO(lll)-0 surface 
fits into this picture, as there the layer distance to the 
topmost oxygen atoms is with 0.51 A significantly smaller 
than for the other orientations, and leads therefore to a 
smaller dipole moment. 

Discarding the polarity as a major factor governing 
the surface stability, we proceed to correlate the latter 
with the binding energy of surface oxygen atoms, as in- 
terestingly each most stable (lxl) termination of each 
orientation features oxygen atoms in the topmost layer, 
cf. Fig. [3 The corresponding binding energies with 
respect to molecular oxygen are shown in Fig. with 
positive values indicating that O2 dissociation would be 



FIG. 5: Binding energies of the topmost oxygen atoms at 
the most stable (1 x 1) termination of the various PdO sur- 
faces. The various orientations are sorted along the x axes 
with higher surface free energies to the right. The binding 
energies are given with respect to a free O2 molecule. 



exothermic. The various orientations are sorted along the 
x axes with higher surface energies to the right, revealing 
a clear correlation. Moreover, there is also a clear correla- 
tion with the coordination of the surface atoms: The two 
most stable terminations, namely the ones contributing 
to our constrained Wulff polyhedron, feature three-fold 
coordinated surface O atoms and with ~ 2.5 eV rather 
strong binding energies. This is followed by the other 
three orientations, that do not contribute to the present 
Wulff shape, exhibiting only two-fold coordinated oxy- 
gen atoms and somewhat lower binding energies around 
1.5 eV. We therefore conclude that the stability of the 
studied subset of (1 x 1) terminations seems primarily 
governed by the openness of the surface orientation, i.e. 
whether its geometric structure offers highly-coordinated 
O binding sites. 



IV. SUMMARY 

In conclusion we have calculated the surface free ener- 
gies of all low-index (lxl) PdO terminations in equilib- 
rium with an oxygen environment using the concept of 
atomistic thermodynamics. The PdO(100)-PdO termina- 
tion exhibits an extraordinarily low surface energy over 
the entire range of experimentally accessible gas phase 
conditions. Correspondingly, this facet dominates the 
Wulff polyhedron constrained to the studied (lxl) ter- 
minations by far, with only the PdO(lOl) orientation 
also covering a smaller surface area. The high stability of 
these two terminations is largely connected to their closed 
geometric structure, allowing a strong oxygen binding in 
highly-coordinated surface sites, while the polarity of the 
non-stoichiometric PdO(100)-PdO termination plays ap- 
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parently only a minor role. The equilibrium shape of a V. ACKNOWLEDGEMENTS 

real PdO crystal is likely to deviate from the presently 
obtained constrained Wulff polyhedron due to surface 
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